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Abstract 

Motivated by recent experiments on the metamagnet FeBr2, anomalies of the 
magnetization and the specific heat in the antiferromagnetic phase of related 
spin models are studied systematically using Monte Carlo simulations. In par- 
ticular, the dependence of the anomalous behavior on competing intralayer 
interactions, the spin value and the Ising-like anisotropy of the Hamiltonian 
is investigated. Results are compared to experimental findings on FeBi"2. 
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I. Introduction 

FeCk* and FeBr 2 are much studied metamagnets of Ising type. The magnetic field 
(if )-temperature (T) phase diagram displays an antiferromagnetically ordered phase, with 
the transition to the paramagnetic phase being of first order at low temperatures and of 
second order at higher temperatures and lower fields. In the antiferromagnetic phase the 
spins of the iron ions are aligned ferromagnetically in the triangular layers perpendicular to 
the c-axis; along that axis there is an antiparallel ordering of the spins. 

In FeC^, the two kinds of transition meet at a tricritical point. For FeBr2, a possible 
decomposition of the tricritical point into a critical end point and a bicritical end point has 
been discussed, in the context of the recent experimental discovery of lines of anomalies in 
the antiferromagnetic phase.iHl3 In particular, the specific heat as well as the temperature 
derivative of the magnetization and the order parameter may display, at fixed field and 
varying temperature, shoulders or maxima below the transition to the paramagnetic phase. 

The anomalies have been attributed&0 to two crucial ingredients of FeBr 2 , the effectively 
weak ferromagnetic intralayer couplings, due to competing antiferromagnetic longer-range 
interactions, and the highly coordinated interlayer couplings to many equivalent iron ions in 
adjacent layers, due to the superexchange mediated by the non-magnetic bromide planes. 
The anomalies have been suggested to reflect the onset of local fluctuations of a second 
antiferromagnetic phase, the All phase, which, if becoming eventually thermally stable, 
would lead to a decomposition of the tricritical point. 

In this article, we shall extend the previous analyses to study quantitatively the depen- 
dence of the anomalies on different parameters of a realistic modeffi for FeBr 2 , namely 
the competing intralayer couplings, the spin value (being 1 in FeBr2), and on the Ising-like 
anisotropies in the Hamiltonian. Furthermore, the relation between the decomposition of 
the tricritical point and the anomalies will be discussed, in particular when comparing our 
Monte Carlo results to recent experimental data and their interpretation. 

The layout of the paper is as follows: The Hamiltonian, obtained from spin wave mea- 
surements, is introduced and experimental findings are outlined in Sect. II. Then results 
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on related Ising models are presented, clarifying the influence of the competing interactions 
and the spin value, followed by a section on the anisotropic Heisenberg model. In Sect. V, 
the comparison to experiments is given. Finally, a brief summary concludes the article. 



II. Realistic Hamiltonian for FeBr 2 

The compound FeBr2 has the hexagonal structure shown in Fig. 1, with the magnetic 
iron ions forming triangular layers perpendicular to the c-axis (corresponding to the z-axis 
of Cartesian coordinates). Based on spin-wave analyses@I!, the low temperature magnetic 
properties of FeBr2 may be obtained from an effective anisotropic Heisenberg Hamiltonian 
for the iron ions 

w = E l-^i^l - Jv { s f s ! + s i s ")) +J2 D {( s ff - ^} - hY.si (i) 

i>j { 'I ) i ^ ' i 

with the spin value S — 1. The first term describes exchange interactions between spins 
in the same triangular layer and adjacent layers. Two different sets of interactions have 
been proposed for the intralayer couplings, with ferromagnetic nearest neighbor interactions, 
J\, and competing antiferromagnetic interactions, extending either up to only next-nearest 
neighbors^, J 2 , 

J l /k B = 7.3 K and J 2 /k B = -2.4 K (2) 
or up to third neighbors in the triangular layersS, J3, 

J x /k B = 4.8 K, J 2 /k B = -0.1 K and J 3 /k B = -1.0 K. (3) 
The interlayer coupling has been determined unambiguously to be 

J[/k B = -2.9 K 

denoting the total exchange to the adjacent iron layer. Taking into account the ten equivalent 
superexchange paths, as mediated by the bromide planes, each individual bond between 
neighboring layers is expected to contribute J'/k B = —0.29 K (see Fig. 2). 



The Ising-type anisotropy, 77 = 0.780, in the first term of the Hamiltonian (1) is enhanced 
by the second term, describing a single-ion anisotropy with the easy axis of the spin along 
the z-axis. Here, D is the energy difference between the doublet and singlet in the lowest 
triplet of an iron ion, with D/ks = —10.7 K (for the intralayer couplings of Eq. (2)) or 
D/ks = —12 K (for the intralayer couplings of Eq. (3)). 

The third term in Eq. (1) describes the effect of the magnetic field, H, applied along the 
c-axis, i.e. in z-direction. 

Fig. 3 shows the H-T phase diagram of FeBr2 determined from measurements of the 
magnetizationi'l, dynamic susceptibility!!, and specific heat. I! Varying temperature, at fixed 
field, all three quantities or their temperature derivatives display in the antiferromagnetic 
phase unusual behavior in the form of shoulders or maxima at about the same temperature 
T a (H), locating the anomaly line. That line seems to evolve from the tricritical point. Note 
that it has been alternately sugg estedB that the anomaly line represents, at sufficiently large 
magnetic fields, a true phase boundary line between different antiferromagnetic orderings. In 
that case, one may expect, from mean-field considerations, the anomaly line to emerge from 
the (bi) critical point at the end of the additional phase boundary line0, with the tricritical 
point having turned into a critical end point. 

In Fig. 3, T p denotes a line in the paramagnetic phase at which the dynamic susceptibility!! 
and the specific heat! show a maximum, when changing temperature at fixed field. It may 
seem to be conceivable that this line also evolves from the tricritical point (or critical end 
point), but this aspect has not been investigated experimentally in detail. 

In the following, we shall study simplified models based on the anisotropic Heisen- 
berg Hamiltonian for FeBr2, Eq. (1), to clarify which of its features may enhance (or 
weaken) the anomalies and, possibly, decompose the tricritical point. So far, previous recent 
analyseJffi dealt with Ising variants of (1), where S = 1/2. Perhaps most importantly, 
the crucial importance of the high interlayer coordination, driving the system close to a 
mean-field type behavior and thereby inducing local thermal excitations of All-type for 
weak intralayer exchange couplings, was established.i'0 Taking merely interactions to the 



geometric nearest neigbor spins in adjacent layers, no anomalies were found&B fin that case, 
fluctuations also destroy the All phase, and hence the tricritical point does not decompose, 
as had been seen in simulation si). Here, we shall elaborate systematically on the role of 
the other parameters in the Hamiltonian, specifically on that of the competing intralayer 
exchange couplings, the spin value, and the spin anisotropy. Thereby, we shall approach 
a rather realistic description of FeBr2. A full analysis of the complete model, Eq. (1), is, 
however, beyond the scope of our study. In addition, such an analysis may provide only an 
integral and thence a rather limited insight into the relevant ingredients leading to anomalies 
in the antiferromagnetic phase of metamagnets. 



III. Ising models 

We shall first approximate the Hamiltonian (1) by Ising models, to elucidate quanti- 
tatively the importance of the competing intralayer couplings as well as the spin value in 
stabilizing the anomalies in the magnetization and the specific heat. 



A. Spin 1/2 

Let us consider the S = 1/2-Ising Hamiltonian 

n = -.r SiSj-j^SiSj-j* £ s^- j 3 £ s^-H^Si (4) 

(NN) NN NNN 3NN i 

where Si is an Ising spin on site i, with spin value 1/2. The exchange interactions Jj.j 
describe, as before, intralayer (extending up to third neighbors in the triangular planes, J\, 
J2, and J3) and interlayer (to the ten equivalent sites in the adjacent plane, J') couplings. 
The couplings are normalized by setting \J'\ — 1. To study the effect of the competing 
interactions in the planes, we usually fix the nearest neighbor interaction J 1; and vary the 
two remaining antiferromagnetic couplings, J2 and J3. According to the two different types 
of exchange constants determined experiment allyBi, two cases are of special interest: (a) 
J 2 = 0, J3 < 0, and (b) J 2 < 0, J3 = 0, respectively. To quantify the efficiency of the antifer- 
romagnetic couplings in weakening the effective ferromagnetic nearest neighbor interactions, 



we also investigated the case (c) J 2 = , J3 = 0, changing J\. 

We simulate systems with K layers, each one consisting of L x L spins, using full periodic 
boundary conditions. Typically, we choose K = L = 20 (to check finite size effects, K and 
L ranged from 10 to 40). For equilibration, 10 4 Monte Carlo steps per site (MCS) were 
used; averages were taken over the following 2 x 10 4 MCS. To improve the statistics and to 
calculate error bars, we performed simulations for ten realizations, with different random 
numbers, at a given field, H, and temperature, T/\J'\. We computed several quantities of 
interest, in particular the energy, E, the specific heat, C (both from energy fluctuations 
and by differentiating the energy with respect to the temperature), the magnetization per 
layer, M(i), and related quantities such as the total magnetization, M, the sublattice mag- 
netizations, Mi and M 2 , referring to the odd and even layers, and the order parameter 
M s = (M\ — M 2 )/2. To take into account phase shifts or flips of entire spin layers, we 
usually computed the absolute values of the total magnetization and the order parameter 
(which will be denoted by M s in the following). In a few cases, we also determined correla- 
tion lengths from standard spin-spin correlation functions. 

In case (a), the ground state, at T = and H < H c q = 20 | J'\, is the antiferromagnetic 
structure, Mi = 1 and M 2 = — 1, assuming |J 3 | < ^Ji (otherwise, more complicated spin 
configurations are stabilizedB, due to the competing interactions along the axes of the tri- 
angular layers). Results of the simulations for that case, fixing the field at H = 0.9 H c0 and 
changing the temperature, are depicted in Fig. 4, showing the specific heat, the order pa- 
rameter, and the temperature derivative of the total magnetization for various values of J3. 
In accordance with the experimental finding for FeBr 2 , J\ has been set equal to 16.75 \ J'\ 
(recall that the values obtained from the spin wave analysis are Jx/ks = 4.8 / 77 = 6.2 K and 
j'/k B = -O.29/77 = -0.37 K). 

In the finite Monte Carlo system, the transition to the paramagnetic phase, at T/v, man- 
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ifests itself, for instance, by a maximum in the specific heat and a drastic decrease in the 
order parameter M s , leading to singularities in the thermodynamic limit. More interestingly, 
anomalous behavior is seen in Fig. 4 to occur well below that transition. For example, the 
specific heat and the temperature derivative of the magnetization display shoulders or max- 
ima, becoming more pronounced with increasing antiferromagmetic interactions J3. The 
anomalies vanish for smaller values of J3 (not shown in Fig. 4). 

Let us briefly recall the physical picture underlying the anomalies, as has been obtained 
from mean-field theory of Ising metamagnets with only nearest neighbor ferromagnetic 
intralayer couplings.0 If those couplings are sufficiently weak, compared to the interlayer 
interactions, a second antiferromagnetic phase, All, may be formed in between the usual anti- 
ferromagnetic phase (AI, with Mi > and M 2 < 0) and the paramagnetic phase (Mi = M 2 ), 
in which both sublattice magnetizations are positive, but different. The All phase may be 
thought of balancing the conflicting tendencies of the external field and the antiferromag- 
netic interlayer couplings, by maintaining, rather small, clusters of 'minus' spins in the even 
layers. Strong ferromagnetic intralayer interactions tend to disfavor those clusters, thereby 
suppressing the All phase. The transition between the AI and All phases is of first order, 
with the boundary line evolving from the critical end point on the border line to the param- 
agnetic phase, and terminating at a (bi)critical point. From that (bi)critical point, a line of 
anomalies emerges. However, such a line may persist even when there is a tri critical point, 
provided the ferromagnetic intralayer couplings are still sufficiently small. 

Including now competing antiferromagnetic intralayer couplings, one may try to cast 
them, together with J 1; in an effective nearest neighbor ferromagnetic interaction, J e ff- To 
elucidate the effect of J 3 on reducing J e ff, we compared our simulational data, case (a), to 
those for models with only nearest-neighbor ferromagnetic intraplane couplings, case (c), 
varying J 1; with \J'\ = 1. In particular, we determined the change in J 1; 5Ji, needed to 
reproduce the Neel temperature T/v, at H = 0, when J3 7^ 0. A naive argument of mean- 
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field type suggests that 8J\ = J3, i.e. J e ff — J\ + J3. In reality, the antiferromagnetic 
coupling is much more efficient in lowering the effective interaction (as may be already seen 
from the analysis of the ground states). For instance, at J\ = 16.75 \ J'\ and J 3 = 4.9 J', 
we find 5Ji ~ —8.15 \ J'\. Note that an even much stronger reduction in Ji is required 
in reproducing, instead of T N , the kind of anomalies present when the antiferromagnetic 
intraplane couplings are included, see Fig. 5. We find, at H = 0.9 H c0 , that the value of 
J 3 = 4.9 J' then corresponds to weakening Ji from 16.75 | J'\ to roughly 1.5 \J'\. The high 
efficiency of J3 (or J 2 , see below) in lowering the effective ferromagnetic coupling J e // and 
hence the ferromagnetic ordering in the layers is crucial, together with the large interlayer 
coordination, in explaining the experimentally found anomalies in FeBr 2 . 

Note that the anomalies shown in Fig. 4 do not correspond to sharp phase transitions. 
For instance, they do not seem to give rise to singularities, as one increases the size of the 
Monte Carlo systems (going from L = K = 10 to 40, e.g., the height of the 'anomalous' 
maximum in C below T N , at H = 0.95 H cQ , does not change significantly, in contrast to the 
behavior of C close to T N) where the peak becomes clearly visible at L = K = 20, increasing 
furthermore for the larger systems), nor is there any indication of hysteresis (by crossing 
the anomalies from different directions in the field-temperature phase plane). Indeed, the 
anomalies may be interpreted as reflecting the onset of local ordering of All-type, with the 
long-range order of the All-phase being, possibly always, destroyed by fluctuations. They 
may be also illustrated by monitoring typical equilibrium Monte Carlo configurations. 

Fixing J 3 and varying the external field, one may map the anomaly line, T a (H). Exam- 
ples for a specific value of J3, J3/J1 = —0.29, chosen to be close, but, in order to identify 
easily the location of the anomalies, somewhat larger than that obtained for FeBr 2 , are 
depicted in Fig. 6. Obviously, the anomalies become stronger upon increasing the field H. 
However, they seem to go over into singularities only at the tricritical point on the phase 
boundary to the paramagnetic phase, as concluded from analyses of the types mentioned 
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above, for fields in the range in between 0.8 H c q and 0.95 H c q (strictly speaking, if there is 
a decomposition of the tricritical point, then the critical end point on the transition line to 
the paramagnetic phase and the critical point at the end of the phase boundary between the 
AI and All phases would be very close to each other). The Monte Carlo data for locating 
the anomaly line are summarized in Fig. 7, depicting the phase diagram in the H-T plane, 
at J3/J1 = —0.29. T a has been determined from the anomaly in the specific heat, in good 
agreement with the corresponding estimates obtained from the magnetizations. 

We also identified the line, T p , in the paramagnetic phase, at which the specific heat C, 
at fixed fields, displays a maximum as a function of temperature, see Fig. 7. The maximum 
is believed to reflect a disordering in the triangular layers.0 In close agreement with the ex- 
perimental findings on FeBr 2 , see Fig. 3, the line seems to evolve from the tricritical point. 
This feature may be, however, accidental. In mean-field theory, T p intersects the boundary 
of the antiferro magnetic phase, T/v, at some point, which is, in general, not related to the 
tricritical or critical end point. In the simulations, the height of the maximum in C does not 
change drastically on approach to the boundary of the antiferromagnetic phase, indicating 
a non-critical behavior, see Fig. 8. It may be wortwhile to clarify this aspect by determining 
the location of T p for different values of J3, where it may be easier to disentangle the inter- 
section point, of T p and T/v, and the tricritical point. Note that the specific heat C, fixing 
the temperature and varying the field, exhibits in the paramagnetic phase a maximum at 
about T p as well, see Fig. 7, in accordance with recent experimental findingsJlI 

In case (b), i.e. J 2 < 0, J3 = 0, similar conclusions hold. To describe experimen- 
tal data on FeBr 2 , we may choose J\ = 25.1 \J'\ and J 2 = 8.4 J'.0 The antiferromag- 
netic couplings lead to a weakening of an effective ferromagnetic intralayer interaction, giv- 
ing eventually rise to All-type excitations in the even or 'minus' layers which cause the 
anomalies in the specific heat C and magnetizations. Actually, J3 is slightly more effi- 
cient than J 2 in reducing J e ff, as seen when adjusting J 2 , with J\ = 25.1 | J'\, to reproduce 
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Tn(H = 0, J\ = 16.75 | J'\, J3 = 4.9 J'). The Neel temperature is realized, when J 2 = 10.3 J' 
(being not far from the experimentally determined value). The anomalies for the two sets 
of parameters do not differ much, see Fig. 5, demonstrating that both types of couplings, J 2 
as well as J3, have a comparable effect on the anomalies, although they are of quite distinct 
physical character (frustration on triangles, J 2 , or competition along the axes of the trian- 
gular planes, J3). It should be emphasized that, in general, frustration or competition is not 
really needed for obtaining the anomalies in the antiferromagnetic phase of metamagnets: 
J e ff has to be sufficiently weak. 

In addition, we determined in which way the ratio of the tricritical temperature T t to the 
Neel temperature Tn(H = 0) depends on the strength of the antiferromagnetic intraplane 
interactions, for the cases (a) and (b), fixing J\ at the value appropriate for FeBr2. The 
ratio decreases with increasing J 2 or J 3 (i.e. decreasing J e ff, see also results from mean-field 
theory, simulations and high-temperature series expansionsi000). For instance, in case 
(a), the ratio varies in between about 0.6 and 0.4, when changing J 3 from 3.3 J' to 6.5 J', 
with Ji = 16.75 \J'\. Similarly, the ratio may be lowered to about 0.32, when increasing J 2 
to 12.1 J', with Ji = 25.1 \ J'\. 

B. Spin 1 

We now consider the S = 1 Ising Hamiltonian, see Eq. (4), where each spin can take the 
values 0, 1 or -1. Compared to the situation with 5* = 1/2, thermal fluctuations are facili- 
tated, reducing the transition temperatures and resulting in more pronounced anomalies. 

In particular, we studied the case Ji = 16.75 \J'\ and J 3 = 4.9 J 7 , with J 2 = D = 0, 
setting I J' I = 1, as before for S = 1/2. Results are displayed in Fig. 9, showing the specific 
heat C versus temperature at various fields, compare to Fig. 6. Clearly, at larger fields a 
maximum shows up below TV; that anomalous behavior in C is corraborated by similar 
properties of the magnetizations, for instance, of dM s /dT. 
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The intersection point of the line of anomalies and the boundary to the paramagnetic 
phase is supposedly the tricritical point (again, we found no evidence for a transition of first 
order between the AI and All phases). The ratio T t /T^{H = 0) is roughly 0.5, as is the case 
for 5 = 1/2 with the same values of J\ and J3. Note that the Neel temperature, Tn(H = 0), 
is, however, compared to its value for the S = 1/2 Ising model, lower by nearly 30 percents. 
Expressing the coupling constants in terms of Kelvin, one easily sees that one moves in 
the case of S = 1 much closer towards the experimentally determined Neel temperature in 
FeBr 2 , see below. 

IV. Anisotropic Heisenberg models 

We now proceed to the anisotropic S = 1 Heisenberg model, given in Eq. (1). In a 
semiclassical description^ of such a model, the z-component of the spin, of lenght 1, is 
discretized, taking the values S z = 0, 1 or -1. If S z = 0, then the spin can rotate, like a 
classical vector, in the xy-plane, see Fig. 10. 

As in the Ising case, we studied especially the case J\ = 16.75 \J'\, J 2 = 0, and J 3 = 4.9 J'. 
Putting D = and rj = 0.78, the simulational data for the specific heat and the magneti- 
zations are very close to those for the corresponding S — 1 Ising model. For example, at 
H = 0.9 H c o, the critical and anomaly temperatures are, in the Heisenberg model, lower by 
roughly one percent. The specific heat is essentially identical to that shown in Fig. 9 for 
the Ising model. In turn, the derivative of the order parameter, dM s /dT, for the Heisenberg 
model, see Fig. 11, agrees very well with that for the Ising case. 

By turning on the single ion anisotropy, D, the critical temperature is shifted towards 
higher values, and the anomalies are somewhat suppressed. In effect, by discriminating 
S z = 0, one approaches the S = 1/2 Ising model. 
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In general, the thermal properties of semiclassical S = 1 Heisenberg models seem to 
resemble quite closely those of the corresponding S = 1 Ising models. Deviations are due to 
spins with vanishing z-component, which may provide, e.g., additional energy contributions. 
Obviously, the different spin components are not decoupled, leading, perhaps, to intriguing 
effects. However, it is beyond the scope of our study, to explore this class of models exten- 
sively (in passing, we may mention our simulational results on the S = 1 Heisenberg model 
with ferromagnetic couplings between neighboring spins on a square lattice. They indicate 
non-critical energy contributions stemming from the xy-components of the spins, leading 
to a minor lowering in the transition temperature, compared to that of the corresponding 
Ising model). 



Note that the discretization of the z-component of the spin is crucial in reproducing 
the anomalous behavior found in FeBr2. A classical Heisenberg model with spins of fixed 
length, but arbitrary orientation, is not expected to show any tendency towards forming, 
even locally, the All phase. Indeed, preliminary simulations on such Heisenberg models did 
not show anomalies in the specific heat or the magnetizations. 



V. Comparison with experiments 

A typical phase diagram of a simplified, but supposedly rather realistic model for FeBr 2 
is depicted in Fig. 7. Obviously, it resembles quite closely the experimental phase diagram, 
see Fig. 3. However, for a quantitative comparison, a few points need to be viewed with 
care. 

Experimentallyffii, the Neel temperature Tn(H = 0) is found to be 14.2 K. Using the 
two sets of coupling parameters as obtained from spin wave analyses, see Eqs. (2) and (3), 
Tn(H = 0) moves towards that temperature from above, by going from the Ising models 
with S = 1/2 to those with S = 1 and finally to the anisotropic Heisenberg models (both sets 
give only slightly different transition temperatures). Indeed, the S = 1/2 Ising models, for 
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both sets of parameters, overestimate Tn(H = 0) by almost a factor of 2 (note that previous 
analyses for FeBr 2 were restricted to that case). For the semiclassical S=l Heisenberg model, 
Tn(H = 0) is about 20 K, i.e. it is still too high. That remaining difference may be partly 
due to a temperature dependence in the effective strength of the single-ion anisotropy, D, 
as had been observed in FeCk^, with D becoming smaller at higher temperatures, thereby 
tending to lower the Neel temperature (in FeBr2, D had been estimated only at a single, low 
temperature).- Similarly, the ratio of T t jT^ is not reproduced quantitatively by the model 
description. While it is about 0.34 in FeBr 2 , the simulations yield such low values, e.g., 
when increasing the antiferromagnetic intraplane interactions beyond the experimentally 
determined values, as discussed above. 

Of course, the deviations from the experimental results might be due to simplifications 
in the model and its treatment, such as neglect of dipolar interactions between the spins 
(their relevance may be seen from the broad two-phase region at low temperatures; they 
also would affect the problem of distinguishing the external, used in experiments, from the 
internal magnetic field, used in the simulations) and neglect of much of the quantum nature 
of the spins. 

As stated before, the main aim of our study is to discuss the origin and character of the 
anomalies in the antiferromagnetic phase. While the model description gives no evidence for 
a sharp transition from the AI to the All phase, such a possibility has been suggested recently 
based on measurements of the specific heat! and, using neutron scattering techniques, the 
order parameter mJI. In particular, the specific heat, as a function of temperature, showed 
a sharp peak superposed on the broad shoulder or maximum well below the transition to the 
paramagnetic phased, becoming sharper with increasing field. However, these findings have 
been questioned lateral. Indeed, no peaks were detected, but only the shoulders or maxima, 
in agreement with the model calculations. 
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In addition, the experimental data for Ms!!! were interpreted in favour of a real transition 
between the AI and All phases. The data, at different fields, were plotted against T/T^(H) 
and against T/T a (if).i In the former case, data separation was observed for T < 0.95 Tjy, 
while in the latter case, the data seemed to fall on one 'universal' curve for T < T a . 

In Fig. 12, we show the corresponding plots of the Monte Carlo data for the S = 1/2 
Ising model with Ji = 16.75 \J'\, J 2 = and J3 = 4.9 J'. Indeed, the behavior is quite 
similar to that found in the experiments, with a clear separation of the order parameter for 
different fields in the predicted ranges of temperatures. Because our analyses, see above, 
give no indication for a sharp phase transition at T a , at least for the fields shown in Fig. 
12, we, however, tend to conclude that this type of data presentation is not suitable for 
providing convincing evidence for the suggested phase transition. 

VI. Summary 

Motivated by recent experiments on the metamagnet FeBr 2 , anomalies in the antiferro- 
magnetic phase of Ising-type models, closely related to the realistic Hamiltonian for that 
magnet as determined from spin wave analyses, have been studied using Monte Carlo tech- 
niques. 

We clarified which ingredients of the Hamiltonian are relevant for the anomalous proper- 
ties, such as broad shoulders or maxima in the specific heat and magnetizations well below 
the transition to the paramagnetic phase. In general, the anomalies can be attributed to 
local thermal excitations of the All phase, due to high coordination of spins in adjacent lay- 
ers and weak effective ferromagnetic intraplane couplings. We demonstrated quantitatively, 
that, in an Ising description, the anomalies are enhanced by competing antiferromagnetic 
intralayer interactions, extending up to third neighbors, and by the spin value, S — 1 in the 
case of FeBr 2 . Going from the Ising model to an anisotropic S = 1 semiclassical Heisenberg 
model leads to only minor changes in the specific heat and the magnetizations parallel to 
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the direction of the applied field. The discretization of the z-component of the spin plays a 
major role in obtaining the anomalies. 

Finally, we compared results of our simulations to experimental findings. The simula- 
tional data suggest that the anomalies usually do not correspond to a sharp phase transition. 
Conflicting interpretations of experiments may be viewed with much care. 
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FIGURES 

FIG. 1. Sketch of the crystal structure of FeBr2, showing the Fe 2+ (full circles) and Br - (open 
circles) ions. 

FIG. 2. The triangular iron planes, with the ten equivalent neighbors (full symbols) in the 
adjacent layer below. 

FIG. 3. Approximate experimental phase diagram, based on measurements of the specific heat 
and magnetizations, see Refs. 5 and 6. T a denotes the anomaly line in the antiferromagnetic 
phase, T p indicates the location of maxima in the specific heat, at fixed magnetic fields, in the 
paramagnetic phase, and Tjy the boundary to the paramagnetic phase. At low temperatures, the 
transition of first order leads to a 2-phase region. 

FIG. 4. Monte Carlo data of (a) the specific heat C, (b) the order parameter M s , and (c) the 
temperature derivative of the total magnetization d\M\/dT versus T/\J'\ for the 5 = 1/2 Ising 
model with J\ = 16.75 \J'\ and H = 0.9 H c q at various values of J3. Systems with K = L = 20 
spins are considered. Here and in the following figures, error bars are only shown when they are 
larger than the sizes of the symbols. 

FIG. 5. Simulational data of the specific heat C versus temperature T/\J'\, at H = 0.9 H c o, 
with (a) Ji = 16.75 |J'|, J 3 = 4.9 J', (b) J x = 25.1 \ J'\, J 2 = 10.3 J', and (c) J x = J eff = 8.6 \ J'\, 
where the Neel temperature Tn(H = 0) is approximately the same in all three cases. Systems with 
K = L = 20 spins are considered. 

FIG. 6. The specific heat C as a function of temperature T/\J'\ for the S = 1/2 Ising model 
with J\ = 16.75 I J' I and J3 = 4.9 \J'\ at various fields. Systems with K = L = 20 spins are 
simulated. 
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FIG. 7. Phase diagram in the field (H)- temperature (T) plane of the S = 1/2 Ising model with 
J\ = 16.75 \J'\ and J3 = 4.9 J'. The anomalies, T a , in the antiferromagnetic phase, determined 
from the specific heat and the magnetizations, are denoted by full circles. The maxima in the 
specific heat ,T p , at fixed fields (open) or fixed temperatures (full), in the paramagnetic phase are 
shown by triangles. Monte Carlo systems with 20 x 20 x 20 spins are simulated. 

FIG. 8. Specific heat C versus temperature, T/\J'\, at various fields in the paramagnetic phase, 
see Fig. 7. 

FIG. 9. Specific heat C versus temperature T/\J'\ for the 5=1 Ising model with K = L = 20 
spins, at 3\ = 16.75 \ J'\ and J3 = 4.9 J', and various fields. 

FIG. 10. Orientations of the spin used in the semiclassical S = 1 Heisenberg model, with 
discretization of the z-component and continuous symmetry in the xy-plane. 

FIG. 11. Temperature derivative of the order parameter, dM s /dT, versus temperature T/\J'\ 
for the anisotropic 5 = 1 Heisenberg model with 20 x 20 x 20 spins, for J\ = 16.75 | J'|, J3 = 4.9 J', 
D = 0, and 7] = 0.78, at H = 0.9 H c0 . 

FIG. 12. Order parameter, M s , versus reduced temperatures, (a) T/Tn(H), and (b) T/T a (H), 
for the S = 1/2 Ising model with J\ = 16.75 \J'\, J3 = 4.9 J', at various fields, showing ranges 
of temperatures (see text) where data (almost) collapse or are widely separated, as in experiments 
on FeBr2, see Ref. 7. Systems with K = L = 20 spins are simulated. 
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Pleimling/Selke: FIGURE 12b 
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